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I.  SUMMARY 
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The  magnetoplasmadynamic  (MPD)  arcjet  has  been  recognized 
over  the  last  decade  and  a  half  as  a  prime  candidate  for  high 
thrust  density,  high  specific  impulse  missions,  such  as  space¬ 
craft  orbital  maneuvering.  Since  the  performance  character¬ 
istics  of  MPD  arc jets  depend  on  proper  matching  of  electro¬ 
magnetic  and  fluid  mechanical  constraints  within  the  thrust 
chamber,  theoretical  modeling  of  the  internal  flow  dynamics 
is  presently  under  development.  Experimental  measurements  of 
internal  flow  conditions  during  MPD  arcjet  operation  are  needed 
to  guide  development  of  flow  models  and  to  assess  the  validity  of 
theoretical  predictions.  Presently,  some  theorists  have  been 
referring  to  data  obtained  by  Turchi  in  1970  applying  rudimentary 
diagnostic  techniques  to  a  quasi-steady  MPD  arcjet  of  non-optimum 
design  (i.e.,  short  cathode,  ablating  insulator,  poorly-mixed 
propellant  feed) .  New  measurements  are  needed  using  modern 
diagnostics  in  MPD  arcjet  flows  that  should  be  more  nearly 
optimal  (based  on  the  last  fifteen  years  of  development) . 

At  the  RDA  Washington  Research  Laboratory,  in  a  cooperative 
effort  with  the  AFOSR-sponsored  MIT  group  under  Prof.  Martinez- 
Sanchez,  an  array  of  diagnostic  techniques  has  been  assembled, 
including  time-,  space-,  and  spectrally-resolved  photography, 
with  which  to  examine  MPD  arcjet  internal  flows.  In  addition, 
the  previous  arcjet  operating  pulsetime  at  RDA  has  been  extended 
from  150  Msec  to  the  millisecond  regime  appropriate  for  studies 
of  megawatt  quasi-steady  and  steady  MPD  arcjets.  The  power 
and  mass  flow  feed  systems  have  been  redesigned  and  rebuilt  to 
accommodate  a  variety  of  MPD  arcjet  geometries  ranging  from 
large  cathode-diameter  annular  channels  with  axial  variation  of 
transverse  dimension  to  channels  formed  by  nested  biconical 
electrodes . 

Efforts  are  continuing  to  apply  the  assembled  apparatus  and 
diagnostic  tools  to  measure  the  electromagnetic  and  plasma  flow 
conditions  in  MPD  arcjets.  Experimental  elucidation  of  the 
internal  flow  structure  will  then  be  used  to  develop  predictive 
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models  for  optimal  geometries  and  operating  parameters.  As  part 
of  the  MIT  modeling  effort,  for  example,  probes  have  been  used  to 
map  the  current  and  voltage  distributions  within  uniform  height 
and  flared  annular  channels.  This  work  is  concerned  with  the 
effect  of  thruster  channel  variations  on  the  current  conduction 
regions  of  the  MPD  internal  flow.  Other  aspects  of  the  internal 
flow  structure  that  have  been  modeled  in  the  past  year  involve 
the  transition  of  the  input  mass  flow  from  neutral  gas  to 
electrically-conducting  plasma.  Scale  sizes  for  the  transition 
region  have  been  estimated  in  terms  of  the  electrical  and 
thermodynamic  properties  of  the  propellant  gas  (e.g.,  argon). 
These  estimates  indicate  that  the  electrical  conductivity  of  the 
flow  can  be  established  in  distances  that  are  small  compared  to 
the  characteristic  dimension  for  current  conduction  near  the 
entrance  to  the  arc jet  thrust  chamber.  Comparison  of  convection 
versus  thermal  conduction  suggests  that  heat  conduction  (and 
associated  surface  ablation)  may  become  more  important  at  higher 
exhaust  speeds.  Examination  of  the  mass  injection  region  will  be 
an  important  area  for  experimental  work  in  the  next  year, 
particularly  in  regard  to  processes  near  the  injectors  that  can 
affect  thruster  lifetime. 


2 


iStk/j 


RDA-TR-144200-001 


SEPTEMBER  1987 


II.  BACKGROUND 

Over  the  past  two  decades,  a  large  number  of  studies  have 
been  performed  in  which  MPD  arc jet  operating  parameters  (e.g., 
current  and  mass  flow  rate)  are  changed  and  variations  in  arc jet 
behavior  are  measured.  The  effects  on  thruster  performance  of 
parameter  variations,  including  overall  electrode  and  injector 
geometry  changes,  have  also  been  noted.  The  basic  problem  is 
that  the  flow  field  and  electromagnetic  structures  within  the  MPD 
thrust  chamber  can  change  significantly  as  operating  conditions 
are  varied,  so  differences  in  performance  may  include  different 
physical  interactions  between  the  internal  flow  and  the  arcjet 
geometry.  Such  differences  in  physical  interaction  include,  for 
example,  enhanced  viscous  effects  along  the  cathode  (center 
conductor)  at  higher  currents  and  lower  mass  flow  rates  simply 
because  the  current  pattern  has  a  greater  axial  component  (due  to 
Hall  effect  and  higher  magnetic  Reynolds  number)  so  the  internal 
flow  is  directed  radially  inward  to  a  greater  degree.  The 
relationship  of  the  mass  injector  positions  to  the  electrode 
surfaces  also  changes  as  variation  of  overall  operating  values 
changes  the  directions  of  current  density  and  acceleration 
patterns . 

The  limited  theoretical  understanding  of  the  MPD  arcjet  and 
its  beginnings  as  a  fortuitously  discovered  mode  of  an  electro¬ 
thermal  arcjet  have  tended  to  freeze  the  electrode  configuration 
in  forms  that  are  unlikely  to  be  ideal.  To  improve  the  MPD 
arcjet  configuration  and  operating  conditions,  however,  guidance 
is  required  in  regard  to  the  effects  of  geometry  and  terminal 
para-  meters  on  internal  processes  and  structures.  Recent 
theoretical  work1,  under  separate  AFOSR  funding,  by  the  MIT  group 
led  by  Martinez-Sanchez  indicates  that  proper  tailoring  of  the 
transverse  dimensions  of  the  thrust  chamber  and  the  distribution 
of  electrode  potential  can  provide  more  uniform  current  distri¬ 
bution  and  possibly  improved  performance  of  MPD  arcjets.  Such 
theoretical  effort  needs  to  be  benchmarked  by  experiments  in 
appropriate  geometries. 
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The  cathode  region  of  an  MPD  arcjet  was  mapped  in  1970  by 
Turchi  using  magnetic  and  electric  field  probes2.  The  electron 
temperatures  were  determined,  along  with  electron  density,  by 
Langmuir  probe  techniques.  On  the  basis  of  this  diagnostic  work, 
Turchi  analyzed  the  electromagnetic  and  flow  structure  of  the 
plasma.  That  this  rudimentary  work,  carried  out  15  years  ago,  is 
still  used  as  a  basic  source  by  people  seeking  to  model  the  MPD 
arcjet  suggests  the  need  for  improved  measurements  with  modern 
diagnostic  techniques.  Furthermore,  the  MPD  thruster  in  which 
this  work  was  done  is,  by  present  standards,  a  non-ideal  configu¬ 
ration.  The  cathode  was  quite  short,  the  gas  flow  ports  were 
poorly  arranged,  and  the  anode  was  essentially  a  ported  planar 
disc.  Useful  results  were  obtained  with  electric  and  magnetic 
measurement  techniques,  but  the  Langmuir  probe  results  for 
temperature  and  density  were  less  reliable.  (Modern  spectroscopic 
diagnostics  are  much  more  appropriate.)  Plasma  flow  characteris¬ 
tics  were  inferred  by  Turchi  using  the  two-fluid  MHD  equations. 
These  measurements  could  be  made  in  a  direct  manner  using  Doppler 
and  spectroscopic  techniques.  Computerized  data  handling  now 
allows  ready  manipulation  of  a  large  number  of  such  measurements, 
so  complete  examination  of  the  MPD  thruster  internal  flow  should 
be  possible. 
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III.  APPROACH 

There  are  two  parts  to  the  present  program  to  develop  a 
useful  understanding  of  MPD  arc jet  internal  flow  dynamics. 
Theoretical  models  are  created  that  predict  relationships  between 
external  arcjet  variables  (e.g.,  geometry,  current,  mass  flow) 
and  internal  structures  (e.g.,  current  conduction  zones,  boundary 
layers,  sheaths) .  Experimental  techniques  can  then  be  applied  to 
examine  the  fluid  mechanical  and  electromagnetic  structure  of  the 
MPD  arcjet  internal  flow  in  order  to  validate  and/or  encourage 
theoretical  modeling. 

Probes  are  still  useful  tools  and  can  be  employed  to  map  the 
current  and  voltage  distributions  in  the  arcjet.  Framing  camera 
photographs  taken  axially  and  through  view  ports  in  the  outer 
electrode  can  be  used  for  survey  work  and  to  identify  flow 
aberrations  (such  as  asymmetries) .  The  electron  density  in  some 
regions  at  least  may  be  deduced  from  spectroscopically-measured 
Stark  broadening  of  selected  spectral  lines.  Temperatures  can 
also  be  measured  using  line  intensity  ratios  and  flow  velocities 
may  be  obtained  from  Doppler  shifts  of  spectral  lines.  In  this 
way,  data  can  be  obtained  on  flow  fields  for  comparison  with  the 
results  of  calculations  using  two- fluid  magnetohydromagnet ic 
theory.  The  experimental  program  and  cross-check  with  theory 
could  be  carried  out  for  a  variety  of  arcjet  configurations 
including  uniform  and  flared  channel  devices  and  for  varying 
conditions  of  current  and  mass  flow  as  appropriate  to  test 
theoretical  understanding  (rather  than  empirically  to  improve 
thruster  performance) .  For  example,  the  fuel  mass  flow  can  be 
reduced  to  provoke  "onset”  conditions  and  the  electromagnetic 
and  flow  characteristics  of  the  internal  flow  associated  with 
this  mode  of  behavior  can  be  examined  directly.  Models  for  the 
MPD  arcjet  internal  flow  are  extended  and  developed  by  using  the 
experimental  results  to  delineate  important  physical  processes, 
such  as  Hall  effect  and  viscous  coupling  to  electrodes  that  might 
otherwise  be  ignored  by  initial  modeling  attempts.  It  is  likely 
that  two-dimensional  MHD  computer  codes  will  be  needed  to 
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calculate  the  basic  flow  accurately,  but  one-dimensional  models 
with  perturbation  techniques  (e.g.,  developed  by  the  MIT  group) 
can  provide  useful  insights  when  combined  with  experimentally- 
determined  regimes  and  modes  of  interest. 
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During  the  first  year  of  the  present  effort,  activities  have 
included  significant  improvement  of  the  facilities  and  diagnos¬ 
tics  needed  for  the  experimental  program,  detailed  examination  of 
the  electromagnetic  structure  in  several  MPD  arc jet  systems,  and 
development  of  models  for  internal  flow  structures.  These 
activities  have  been  a  cooperative  effort  by  the  RDA  Washington 
Research  Laboratory  and  the  MIT  group  under  Prof.  M.  Martinez- 
Sanchez,  under  separate  AFOSR  sponsorship. 


Experimental  Efforts 


In  this  cooperative  effort,  the  power  supply  was  upgraded 
from  a  20  kJ,  150  /isec  pulse  forming  network  to  a  400  kJ, 

0.5-1  msec  LC-ladder  style  pulseline.  Appropriate  changes  in 
the  current  transmission  lines  and  coupling  to  the  arc jet  chamber 
and  in  the  arc jet  electrical  feed  were  also  accomplished.  The 
arcjet  feed  design  allows  several  types  of  arcjets  to  be  tested 
without  changing  the  basic  power  and  gas  handling  connections. 

In  addition,  the  gas  feed  was  upgraded  to  match  better  to  the 
higher  currents  (50  kA  vs  12  kA)  provided  by  the  new  power 
source. 


In  addition  to  probe  techniques,  optical  diagnostics  have 
been  developed  and  improved.  Spatial  distributions  of  electron 
density  have  been  determined  by  measuring  the  Stark  broadening 
profiles  of  the  observed  emission  lines.  With  the  same  spectro¬ 
scopic  system,  plasma  flow  velocities  will  be  obtained  from 
measurement  of  the  Doppler  shift  of  spectral  lines.  The  spectral 
analysis  system  now  available  can  measure  a  shift  of  0.006  nm 
using  third  order  Ar  II  line,  which  implies  a  minimum  detectable 
velocity  of  approximately  3700  m/s.  Plasma  electron  temperature 
will  be  determined  by  the  Boltzmann  plot  technique  in  which 
relative  intensities  of  spectral  lines  are  compared. 

To  enable  spectroscopic  analysis  of  the  entire  flow  field,  a 
digital  spectral  and  spatial  data  acquisition  system  has  been 
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developed  comprising  a  1.2  meter,  f/11.5  spectrograph  with  a 
EG&G/PARC  intensified  optical  multichannel  analyzer  (SIT 
camera) .  The  spectral  and  spatial  output  from  the  camera  is 
processed  by  an  IBM  System  9000  laboratory  computer  for  immediate  3 

background  subtraction  and  signal  averaging.  Extensive  computer 
software  is  available  to  smooth  the  data,  perform  background 
subtraction,  and  line  intensity  integration.  An  Abel  inversion 
routine  can  then  be  used  to  yield  line  profile  measurements  or 
spectral  line  intensities  for  the  axisymmetric  plasma  flow.  A 
more  detailed  description  of  the  digital  spectral  and  spatial 
data  acquisition  system  is  provided  in  Appendix  I. 

Large  diameter,  annular  arc jets  have  been  attempted  in  order 
to  achieve  better  comparison  with  theoretical  models  that  do  not 
include  cylindrical  effects.  The  length  of  the  channel  has  been 
adjusted  and  the  total  current  increased  to  obtain  acceptable 
discharge  uniformity.  Flared  channels  have  been  compared  with 
channels  of  uniform  height  (interelectrode  gap  vs  axial  position) 
to  check  theoretical  predictions  of  current  density  distribu¬ 
tions.  Evidence  of  substantial  axial  current  in  the  flow  (Hall 
effect)  indicates  the  need  for  additional  modeling  activity. 

Preliminary  results  of  experiments  with  large  diameter  annular 
arc jets  were  presented  at  the  19th  International  Electric 
Propulsion  Conference  (Appendix  II) .  More  detailed  examination 
will  be  the  subject  of  a  doctoral  dissertation  by  D.J.  Heimer- 
dinger  (to  be  reported  separately  under  AFOSR-sponsorship) . 

Theoretical  Efforts 

In  a  self-field  magnetoplasmadynamic  arcjet,  the  main 
structure  of  the  discharge  is  governed  by  the  competition  between 
convection  and  diffusion  of  magnetic  flux.  For  high  magnetic 
Reynolds  number,  Rm  =  anu^L,  this  competition  favors  convection 
and  can  result  in  a  bifurcated  current  distribution.  The  current 
density  is  highest  in  regions  near  the  entrance  and  exit  of  the 
MPD  channel,  where  the  back  electromotive  forces  (EMF)  are 
lowest,  (at  the  entrance  because  the  flow  speed  is  low,  and  at 
the  exit  because  the  magnetic  field  must  decrease  to  zero) . 
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Analysis1  based  on  high  magnetic  Reynolds  number  flow,  Rm  »  1, 
suggests  that  the  axial  distribution  of  current  can  be  controlled 
by  varying  the  interelectrode  gap  along  the  channel .  Other  work2 
utilizes  the  magnetic  Reynolds  number  to  scale  the  axial  thick¬ 
ness  of  the  main  current  conduction  regions  of  the  channel, 
thereby  obtaining  a  relationship  between  heat  generation  (due  to 
resistive  dissipation)  and  thrust  power.  This  relationship  is 
then  used  to  examine  the  behavior  of  high  power  MPD  arc jets  in 
terms  of  voltage-current  characteristics  and  possible  exhaust 
velocity  limitations  (Appendix  III)  . 

It  is  assumed  in  the  preceding  analyses  that  the  plasma  flow 
has  achieved  a  sufficient  (and  uniform)  electrical  conductivity 
to  justify  the  high  Rm  approximation.  The  initial  temperature  of 
the  propellant  mass  flow,  however,  as  it  enters  the  MPD  thruster 
is  low  (=  room  temperature)  so  its  electrical  conductivity  is 
negligible.  It  is  useful  to  consider  the  transition  of  the 
propellant  flow  from  very  low  to  high  conductivity  due  to  resis¬ 
tive  heating.  Since  the  geometry  of  the  injector  region  may 
place  this  transition  region  electrically  in  parallel  with  the 
MPD  discharge  flow,  a  connection  might  be  expected  between  the 
propellant  electrical  and  thermodynamic  properties  and  the 
overall  power  distribution  in  the  arcjet.  A  sub-structure  of 
the  discharge  flow  at  the  MPD  channel  entrance  can  be  delineated, 
perhaps  resembling  the  flame  zone  in  propellant  combustion. 

The  size  of  this  structure  may  be  compared  with  the  dimensions 
for  other  processes  to  suggest  possible  consequences  for  MPD 
arcjet  operation. 

Current  Integral  Approach 

In  the  analysis  of  electrically-exploded  metallic  fuses, 
such  as  wires  or  foils,  it  is  common  to  invoke  the  change  in 
electrical  properties  with  change  in  thermodynamic  state  of  the 
fuse  material.3  Typically,  the  operation  of  high  power  fuses 
occurs  on  a  time  short  compared  to  that  required  for  heat 
transfer  (as  opposed  to  lower  power  fuses  that  follow  Preece's 
law4).  The  temperature  and  state  of  the  fuse  thus  depends  on  the 
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heat  deposited  by  resistive  dissipation: 


=pcdT 


(1) 


where  a  =  electrical  conductivity,  j  is  the  instantaneous 
current  density,  p  *  mass  density,  c  =  heat  capacity  per  unit 
mass,  and  the  elapsed  time  and  temperature  increase  are  respec¬ 
tively,  t  and  T.  For  many  materials,  under  a  range  of  condi¬ 
tions,  the  electrical  conductivity  will  be  a  function  of 
temperature.  Equation  (1)  can  then  be  rearranged  to  form  the 
so-called  current  integral: 


J  j2dt  =  |  a p cdT  ,  (2) 

where  the  right-hand  side  depends  only  on  material  proper¬ 
ties  and  may  be  integrated  between  particular  states,  such  as 
room  temperature  to  vaporization.  The  left-hand  side  does  not 
depend  on  the  instantaneous  circuit  current,  J,  but  may  be 
adjusted  by  choosing  the  cross-sectional  area  of  the  fuse  to 
achieve  a  current  density,  j,  such  that  integration  over  the 
circuit  waveform  provides  the  desired  fusing  action  after  a 
specified  time.  Fuses  designed  in  this  manner  have  been  quite 
successful  in  a  variety  of  experiments.5 

Near  the  entrance  to  the  MPD  channel,  where  the  flow 
velocity  is  much  lower  than  the  eventual  exhaust  speed,  it  is 
possible  to  apply  the  current  integral  approach  to  convecting 
fluid  elements  and  thereby  estimate  the  distance  required  to 
achieve  high  electrical  conductivity.  If  heat  transfer  and 
conversion  of  thermal  energy  to  kinetic  energy  are  neglected 
as  initial  approximations,  then  a  fluid  element  of  fixed  mass 
will  experience  a  temperature  increase  due  to  resistive  heating 
as  it  enters  the  main  discharge  region: 

pCpdT  =  <7E2dt  (3) 
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where  E  is  the  electric  field  across  the  main  discharge. 
This  electric  field,  in  a  steady  or  quasi-steady  arcjet,  must 
also  be  supported  by  the  propellant  flow  (downstream  of  the 
injector  plate)  as  it  approaches  and  enters  the  main  discharge. 

The  time  spent  by  the  fluid  element  at  any  position  x  to 
x  +  dx  is  merely  dt  =  dx/u,  where  u  is  the  flow  speed  at  that 
point.  The  local  flow  speed,  however,  may  be  written  in  terms 
of  the  mass  flow  per  unit  area  ma  and  the  local  mass  density 
u  =  ma/p ,  so 


(4) 


For  a  one-dimensional  flow  ma  is  independent  of  x,  so  integration 
of  equation  (4)  provides  a  distance  a  proportional  to  an  integral 
of  flow  properties  over  the  temperature  interval: 

ma  r  cp 

A  =  I?  J  T  dT  *  (5) 

This  calculation  is  completely  analogous  to  that  for  the  explo¬ 
ding  fuse,  with  the  caution  that  the  conductivity  and  specific 
heat  of  a  partially  ionized  gas  are  not  functions  simply  of 
temperature;  in  the  range  of  0.5  to  2  eV,  however,  calculated 
properties  of  argon,6'7  for  example,  suggest  that  the  ratio  of 
enthalpy  to  conductivity  only  varies  by  about  a  factor  of  two 
for  a  pressure  variation  of  two  orders  of  magnitude  (0.01  to 
1.  atmospheres) .  The  present  discussion  is  primarily  concerned 
with  scaling  of  discharge  structures,  which  should  not  be 
obscured  by  absolute  values  that  will  vary  with  the  detailed  flow 
chemistry. 

The  resistive  heat  integral 


H  (T) 


(6) 


is  displayed  in  Figure  1  as  the  scaled  thickness  E2A/ma  =  H(T)  vs 
the  electrical  conductivity  a  .  An  order  of  magnitude  decrease  in 
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conductivity  occurs  from  T  »  20,000*  K  down  to  T  -  6,000*  K,  and 
thereby  defines  the  effective  thickness  of  the  transition  region 
between  non-conducting  and  highly  conducting  flow.  With 
numerical  values  typical  of  experimental  operation  (E  *  3500  V/m 
and  ma  =  0.5  kg/m2-s) ,  the  actual  transition  region  thickness 
for  argon  would  be  about  1  mm,  which  is  much  less  than  tlu 
characteristic  discharge  thickness  based  on  the  magnetic  Reynolds 
number  (5  =  l/onu^  =1-2  cm) . 

Scaling  of  the  Transition  Region 

For  a  given  value  Hf  =  H(Tf),  corresponding  to  the  full 
operating  conductivity  of  the  MPD  flow,  the  size  of  the  transi¬ 
tion  region  will  depend  on  the  electric  field  established  across 
the  main  discharge.  This  electric  field  is  directly  related  to 
the  back  EMF  in  the  high  magnetic  Reynolds  number  flow,  which 
scales  as  the  product  of  exhaust  speed  U£  and  magnetic  field  B. 
Thus: 

E  —  fU£B  ^  7  j 

where  the  value  of  f  depends  on  the  geometry  of  the 
thruster;  (since  ug  and  B  are  the  maximum  values  of  these 
quantities  in  the  flow,  f  will  be  less  than  unity,  but  must 
also  be  adjusted  to  correspond  to  the  appropriate  electric  field 
at  the  radius  of  the  mass  flow  injection) .  The  transition 
thickness  is  then: 

maHf 

A  =  — — ~  r  .  (8) 

f2U£2B2 

The  maximum  magnetic  field  is  proportional  to  the  total 


current 


and  the  exhaust  velocity  for  the  MPD  thruster,  based  on  electro¬ 
magnetic  forces  may  be  written  as 
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mJ2  g 

ue  "  aT  £  <1( 

where  g  is  a  factor  depending  on  the  geometry  of  the 

thruster,  for  example  g  »  (In  rpyrc  +  3/4) ,  based  on  the  ratio 

of  anode  and  cathode  radii  r^/rc.  Note  that  the  full  mass  flow 
•  • 

m  -  maA  is  utilized  here.  Substitution  and  rearrangement  of 
terms  in  equation  (8)  then  provides: 

G  Hf 


The  factor  G  includes  all  the  previous  terms  related  to  the 
thruster  geometry  (and  also  other  physical  constants) . 

Equation  (11)  indicates  that  the  transition  thickness 
decreases  rapidly  with  increased  exhaust  speed.  In  fact,  since 
the  main  discharge  thickness  scales  as  u^"1,  the  transition 
region  becomes  a  progressively  smaller  portion  of  the  discharge 
as  J2/m  increases.  With  the  earlier  numerical  estimate  indica¬ 
ting  a  transition  thickness  that  is  less  than  ten  percent  of  the 
magnetic  Reynolds  number  thickness,  it  appears  that  approximating 
the  MPD  flow  with  a  uniform,  high  electrical  conductivity  is 
indeed  valid  for  the  entrance  region.  Even  though  the  main 
discharge  thickness  decreases  with  higher  magnetic  Reynolds 
number,  the  establishment  of  a  high  conductivity  flow  occurs 
on  still  smaller  dimensions.  A  limit  to  this  benign  situation 
occurs,  however,  due  to  heat  conduction. 

Thermal  Conduction  in  the  Transition  Region 

As  the  dimension  of  the  transition  region  decreases,  the 
temperature  gradient  in  the  flow  increases.  At  some  point,  the 
initial  model  in  which  resistive  dissipation  simply  increased 
the  flow  temperature  becomes  totally  inadequate.  It  is  useful 
to  examine  the  opposite  limit,  namely,  a  balance  between  heat 
conduction  and  resistive  dissipation: 


dr  3T^  - 

-  (k(T)  +  o  (T)  E2  =  0. 
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Both  the  thermal  conductivity  K(T)  and  the  electrical  conduc¬ 
tivity  are  rapidly  varying  functions  of  temperature  in  the 
transition  region.  Inspection  of  the  behavior  of  thermal  and 
electrical  conductivity  as  functions  of  temperature  and  pressure, 
calculated  in  Ref  7,  suggests  that  in  the  transition  region  the 
following  forms  may  be  adequate  for  the  present  discussion. 

K  =  K0e“Ck/T  and  a  -  <70e'Cc/T  .  (13) 

Such  formulation  allows  the  simplification: 


where  w  =  C^Cj,  reflecting  the  typical  situation  that  good 
electrical  conductors  are  good  thermal  conductors.  Equation  (12) 
may  then  be  nondimensionalized  in  terms  of  normalized  variables 

c»  =  X/XQ  ,  S  *  a/aQ  (15) 

so  that 


a 

sw-1  as’ 

fa0E2X02'| 

da 

(In  S)2  da 

KqCq 

With  F(s)  =  sw"V(ln  s)2,  this  equation  simplifies  to: 


(16) 


(17) 


when  the  characteristic  distance  for  change  in  electrical 
conductivity  is  identified  as: 


x0 


K0Cc 

-  a  o  - 


/E 


(18) 


In  Figure  2,  a  solution  of  Equation  (17)  is  displayed  for  the 
case  of  s  =  0.1  in  the  main  discharge.  Note  that  an  order  of 
magnitude  change  in  electrical  conductivity  requires  a  distance 
of  about  2x0.  To  match  equation  (13)  to  the  results  in  Ref  7, 
appropriate  values  for  the  transition  region  in  an  MPD  arc jet 
using  argon  would  be  K0  =  5  w/m-*K,  =  1.5  x  105  mho/m,  and 
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Cc  =  4.6  x  104  *K.  For  E  =  3500  V/m,  as  before,  the  charac¬ 
teristic  scale  size  below  which  thermal  conduction  will  tend 
to  prevent  resistive  heating  of  the  input  mass  flow  to  high 
conductivity  is  then  xQ  =  0.3  mm.  At  the  specified  value  of 
electric  field,  xQ  <  a ,  but  a  factor  of  two  increase  in  electric 
field  (at  fixed  ma)  could  reverse  this  inequality. 

In  summary,  it  appears  that  the  size  of  the  transition 
region  is  much  smaller  than  typical  MPD  discharge  dimensions,  but 
could  become  comparable  to  gradient  scales  for  heat  conduction  as 
the  electric  field  supported  by  the  discharge  increases  with 
higher  exhaust  velocity.  Such  higher  velocity  operation  reduces 
the  thickness  of  the  MPD  discharge,  increasing  resistive  dissipa¬ 
tion.2  Perhaps  more  importantly,  it  allows  the  high  temperature, 
high  conductivity  working  fluid  to  be  generated  closer  to  the 
injection  region.  This  proximity  increases  the  heat  conduction 
to  cold  boundaries  removing  energy  that  would  otherwise  be  used 
more  efficiently  to  create  the  plasma  flow.  In  some  instances, 
the  increased  heat  flux  will  cause  surface  ablation,  resulting 
in  longterm  thruster  degradation  and  lower  exhaust  speeds  (as 
discussed  in  Ref  2) .  Higher  power  operation  can  thus  result 
in  less  performance  than  anticipated  because  of  increased  heat 
transfer  and  its  consequences.  Since  the  transition  region 
electrically  shunts  the  MPD  discharge,  more  detailed  study  of 
the  injector  region  is  probably  warranted  in  order  to  understand 
the  effects  of  injector  and  backplate  design  on  MPD  arcjet 
behavior. 
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V.  CONCLUDING  Ri 


Consideration  of  the  MPD  arcjet  internal  flow  from  a  quasi- 
one  dimensional  standpoint  provides  structural  details  that 
manifest  the  competition  among  processes  of  convection  and 
diffusion.  The  current  density  distribution  for  an  acceleration 
channel  of  uniform  transverse  dimension  will  tend  to  bifurcate 
into  high  current  density  regions  at  the  entrance  and  exit  of  the 
channel.  The  latter  concentration  in  a  one-dimensional  formula¬ 
tion  depends  on  the  requirement  of  zero  magnetic  field  at  the  end 
of  the  channel.  This  requirement  is  artificial  for  a  real  MPD 
accelerator  exhausting  into  vacuum,  and  may  be  replaced  by  two- 
dimensional  flow  regions  extending  from  the  edges  of  the  channel 
exit.  In  the  high  magnetic  Reynolds  number  limit  (at  high 
density) ,  these  regions  (for  rectangular  geometry)  are  centered- 
expansion  fans,  with  characteristic  lines  corresponding  to 
current  streamlines.  Current  concentrations  thus  occur  at  both 
electrode  edges  (for  a  rectangular,  parallel-plate  accelerator) . 
In  coaxial  geometry,  there  is  a  current  concentration  along  the 
inside  diameter  of  the  anode  orifice,  and  also  downstream  of  the 
end  of  the  cathode  (where  the  radially-inward  expansion  of  the 
annular  MPD  flow  corresponds  to  a  non-simple  region  of  interact¬ 
ing  expansion  fans  and  reflected  compression  waves) . 

The  current  concentration  at  the  entrance  to  the  MPD 
thruster  occurs  where  the  electrically-conducting  flow  has  a 
relatively  low  speed  (and  thus  acts  as  a  low- impedance  shunt 
across  the  linear  electric  motor  represented  by  the  thruster) . 
Magnetic  field  can  diffuse  into  this  low  speed  flow  before 
subsequent  acceleration  increases  the  flow  speed  sufficiently 
to  convect  magnetic  flux  downstream.  A  scaling  analysis 
(Appendix  III)  indicates  that  higher  flow  speed  increases  the 
current  concentration  at  the  entrance  resulting  in  greater 
resistive  dissipation  and  increased  heat  transfer  to  the  arcjet 
electrodes  and  back  plate.  As  a  first  step  in  examining  the 
energy  transport  near  the  arcjet  entrance,  the  transition  of  the 
input  mass  flow  from  neutral  gas  to  electrically-conducting 
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plasma  was  modeled  during  the  past  year.  Further  work  needs 
to  be  done  to  consider  the  mass  injection  region  in  two  (and 
possibly  three)  dimensions  in  order  to  understand  properly 
the  upstream  boundary  conditions  that  this  region  may  provide 
in  overall  MPD  arcjet  operation,  performance  and  lifetime. 
Detailed  experimental  examination  of  the  injector  region  will 
be  very  important  in  developing  appropriate  models.  (Older 
works  at  Princeton  by  Turchi  with  probes  and  spectroscopically 
by  A.  Bruckner9  indicate  that  near  the  injectors  there  are 
several  interesting  features  that  can  be  successfully  diagnosed.) 

Since  the  overall  flow  field  within  the  MPD  arcjet  thrust 
chamber  is  not  quasi-one  dimensional,  the  placement  of  the  mass 
injectors  relative  to  the  current  distribution  can  change  as 
current  and  mass  flow  values  are  varied.  Thus,  the  effective 
geometry  of  the  device  changes  with  operating  conditions  and 
must  be  measured  (spectroscopically  and  with  probes)  in  order 
to  understand  arcjet  behavior.  The  MPD  flow  regime  may  be 
inherently  two-dimensional,  so  the  need  to  establish  by  experi¬ 
mental  diagnosis  the  flow  structure  for  each  set  of  operating 
conditions  is  probably  unavoidable. 

The  two-dimensionality  of  the  MPD  flow  field  further  compli¬ 
cates  comparison  of  experimental  and  theoretical  results  because 
simple  solutions  appropriate  to  the  basic  internal  flow  do  not 
exist.  (For  example,  in  the  MPD  regime,  there  is  a  significant 
axial  component  to  the  current,  so  the  main  flow  is  accelerated 
to  high  Mach  numbers  at  finite  angle  of  attack  onto  the  cylindri¬ 
cal  blunt  body  represented  by  the  arcjet  cathode.)  Development 
and  application  of  two-dimensional  MHD  code  techniques  will  be 
necessary  in  order  to  progress  beyond  qualitative  guidance  in  the 
design  of  MPD  arcjets.  Some  work  has  started  in  this  direction 
at  MIT,  under  AFOSR-sponsorship.  On  other  research  programs, 
RDA/WRL  has  been  closely  involved  with  two-dimensional  MHD  codes 
for  coaxial  plasma  guns  (AFWL)  and  dense  plasmajets  (AFOSR/SDIO) , 
so  there  may  be  opportunities  to  apply  existing  capabilities  to 
the  MPD  arcjet  problem. 
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APPENDIX  I 

Details  of  Spectroscopic  Analysis  System 

To  enable  spectroscopic  analysis  of  the  entire  flow  field, 
a  digital  spectral  and  spatial  data  acquisition  system  has  been 
developed  in  a  cooperative  effort  with  the  AFOSR-sponsored  MIT 
group  under  Prof.  Martinez-Sanchez .  The  system  is  comprised  of 
a  1.2  meter,  f/11.5  spectrograph  with  an  EG&G/PARC  intensified 
optical  multichannel  analyzer  (SIT  camera) .  The  spectral  and 
spatial  output  from  the  camera  is  processed  by  an  IBM  System  9000 
laboratory  computer  for  immediate  background  subtraction  and 
signal  averaging.  The  digital  spectral  and  spatial  data 
acquisition  system  development  consists  of  interfacing  systems, 
characterization  of  SIT  camera  performance,  relative  spectral 
calibration,  optical  alignment  with  spatial  calibration, 
development  of  gate  operation  capability  and  computer  analysis  of 
data  with  associated  software  development.  Also,  considerable 
effort  was  necessary  to  shield  the  SIT  camera  and  associated 
electronics  from  electromagnetic  noise  generated  by  the  arcjet 
power  supply  and  connections.  Many  issues  as  to  SIT  system 
performance  surfaced  during  system  characterization  and  calibra¬ 
tion.  These  include  linearity  of  signal,  spatial  uniformity  of 
responsivity ,  background  noise  due  to  internal  electronics,  shot 
noise  due  to  spurious  electrons  from  the  photocathode,  and  edge 
of  field  effects  that  cause  pincushion  distortion  at  frame  edges 
(straight  lines  are  observed  as  curves) .  More  detailed  discus¬ 
sion  of  the  digital  spectral  and  spatial  data  acquisition  system 
and  analysis  of  data  from  large  diameter,  annular  arcjets  will  be 
the  subject  of  a  master's  thesis  by  Daniel  Kilfoyle  (to  be 
reported  separately  under  AFOSR-sponsorship) . 

The  SIT  camera  has  a  detector  array  of  512  x  512  pixels 
(262k  pixels)  with  spatial  size  of  25  microns  and  a  dynamic 
range  of  16k.  The  interfaced  system  had  resolution  less  than 
expected  due  to  signal  blooming  from  pixel  to  pixel  (thus 
decreasing  the  observed  spatial  resolution)  and  range  less  than 
expected  due  to  camera  electronics  (e.g.,  A/D  converters)  and  to 
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software  capabilities.  The  detector  system  was  therefore  limited 
to  250  horizontal  pixels  (spectral  data)  and  to  128  vertical 
pixels  (spatial  data)  or  32k  data  points  per  frame.  The  detector 
showed  excellent  spatial  uniformity  of  responsivity  but  limited 
range  of  signal  linearity  of  0  to  600  counts  as  shown  in  Figure 
1-1.  The  detector  system  relative  spectral  response  was  measured 
with  a  calibrated  tungsten  filament  lamp.  The  lamp  was  mounted 
to  include  effects  of  the  imput  optics,  turning  mirrors  and 
attenuation  due  to  the  chamber  window  material.  The  resultant 
relative  spectral  response  normalized  to  1.0  at  700  nm  is  shown 
in  Figure  1-2. 

For  initial  system  characterization,  collimating  input 
optics  with  a  magnification  of  5:1  (5  cm  at  source)  and  2  cm 
downstream  of  an  arc jet  were  used.  Figure  1-3  shows  raw  data  as 
collected  by  the  SIT  camera  system  with  the  centerline  of  an 
argon  arcjet  at  row  38  (approximately  1/3  up  from  the  bottom  of 
the  page).  Each  vertical  pixel  represents  approximately  0.05  cm. 
This  raw  data  set  has  been  zero-corrected,  which  effectively 
removes  background  noise  due  to  the  SIT  internal  electronics. 

Shot  noise  that  is  random  in  time  for  each  pixel  produces  a 
random  signal  fluctuation  on  the  detector  spatially  but  with 
spatial  frequencies  on  the  order  of  0.33  to  0.5  cycles  per  pixel. 
Figure  1-4  shows  a  hortizontal  scan  of  Figure  1-3  of  row  55  (ap¬ 
proximately  0.85  cm  from  the  experimental  center  line).  Note 
that  each  frame  can  contain  spatial  intensity  information  of  more 
than  one  plasma  emission  line,  in  this  case  emission  from  Ar  II 
lines  460.956  nm,  458.90  nm  and  457.935  nm. 

Software  has  been  developed  for  further  data  reduction  and 
analysis.  For  spatial  electron  density  measurements  the  line 
emission  profiles  is  measured  and  for  spatial  electron  tempera¬ 
ture  determination  the  relative  line  emissivity  is  measured. 

Plasma  electron  temperature  will  be  determined  by  the  Boltzmann 
plot  technique  in  which  relative  intensities  of  several  spectral 
lines  are  compared.  Figure  1-5  shows  row  32  of  an  Ar  I  emission 
line  that  has  been  zero-corrected.  A  small  buffer  containing 
only  the  line  spatial  emission  is  made  and  integrated  to  yield  a 
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Intensity  Response  Curve 


Actual  Intensity 

Figure  1-1.  SIT  camera  observed  counts  (vertical  axis)  vs 

relative  input  light  flux  (horizontal  axis)  at  two 
horizontal  rows  shows  good  spatial  uniformity  of 
responsivity  and  linear  response  to  600  counts. 


Wsvalangth  (no) 

Figure  1-2.  Relative  spectral  system  response  normalized  to  1.0 
at  700  nm  as  measured  by  calibrated  tungsten 
filament  lamp  between  400  and  800  nm. 


Figure  1-3.  Visual  representation  of  SIT  camera  data  of  an  argon  arc jet. 

Each  SIT  camera  data  frame  contains  information  in  three 
dimensions:  spectral,  spatial,  and  intensity. 
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raw  intensity  spatial  profile  shown  in  Figure  I-6a.  This  data  is 
then  digitially-smoothed  via  a  fast  Fourier  transform  technique 
to  eliminate  higher  spatial  frequencies  (shot  noise)  with  result 
as  shown  in  Figure  I-6b.  This  smoothed  raw  intensity  spatial 
profile  is  then  Abel-inverted  to  yield  spatially-resolved, 
spectral  line  emissivity  that  can  be  used  for  input  to  a 
Boltzmann  plot  routine  to  determine  spatially-resolved  electron 
temperature,  as  shown  in  Figure  1-7. 

Figures  1-8  and  1-9  show  zero  corrected  data  for  656.285  nm 
hydrogen  emission  2  cm  downstream  of  the  arcjet  with  1.6% 
hydrogen  added  to  the  argon  in  the  plenum.  For  determination  of 
electron  density,  the  emission  line  profile  and  FWHM  (full  width 
at  half  maximum)  must  be  measured.  Software  manipulation  of 
emission  line  profile  data  is  therefore  performed  both  spatially 
and  spectrally.  Figures  I-10a  and  I-10b  show  the  hydrogen 
emission  line  at  one  radial  position  before  and  after  smoothing. 
Figures  I-lla  and  I-llb  show  unsmoothed  and  smoothed  spatial 
emission  profile  along  the  hydrogen  emission  peak.  The  smoothed 
data  is  then  Abel-inverted  to  yield  true  local  line  profiles. 

This  information  is  first  corrected  by  subtraction  of  line 
broadening  due  to  temperature  (Doppler  broadening)  and  due  to 
instrument  broadening.  The  data  is  then  fitted  to  a  Lorenztian 
profile  to  yield  spatialy-resolved  electron  density. 

Full  flow-field  visualization  of  many  of  the  important  para¬ 
meters  has  been  demonstrated  with  the  new  digital  data  acquisi¬ 
tion  system  allowing  more  detailed  and  accurate  data  to  be 
obtained  in  less  time.  Future  system  improvements  are  expected 
to  increase  sensitivity  and  dynamic  range.  In  addition,  software 
is  being  constantly  upgraded  for  easier  data  reduction  and 
presentation. 
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Figure  1-7.  Abel  inversion  of  digitally-smoothed  integrated  argon  emission 
shown  in  Figures  1-5,  I-6a  and  I-6b  yield  spatially-resolved 
spectral  line  emissivity. 


Figure  1-8.  Visual  representation  of  656.285  nm  hydrogen  emission  2  cm 
downstream  of  arcjet  exit  as  observed  by  SIT  camera  system. 
Analysis  of  hydrogen  spectral  profile  yields  spatial  electron 
density . 


Top:  1380  Basel inn:  1  »Y=  42  with  250  points. 


igure  1-9.  Intensity  vs  wavelength  of  656.285  nm  hydrogen 
downstream  of  arcjet  exit  at  one  spatial  locat 


pectral  emission  line  profile  of  656.285  nra  hydrogen  emission 
rom  data  at  2  cm  downstream  of  an  argon  arcjet  exit  plane 
1.6%  hydrogen  added  to  argon).  (a)  Raw  data  emission  profile; 
b)  digitally-smoothed  emission  profile. 
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Abstract 

Two  annular  magnetoplasmadynamic  MPD 

thrusters  were  constructed  and  operated  in  a 
regime  where  anode  starvation  processes  limit  the 
thrusters'  performance.  Magnetic  field  and 

floating  potential  measurements  were  conducted  to 
determine  the  operational  differences  between  a 
-nannel  with  constant  mtereiectrode  separation 
and  a  channel  with  an  initial  constant  separation 
and  a  iivergent  exit.  The  channel  with  the 
onstant  electrode  separation  shows  evidence  of 
strong  concentrations  of  current  at  both  inlet  and 
exit,  along  with  strong  indications  of 
accompanying  electrode  erosion.  Tn  the  divergent 
section,  current  densities  were  greatly  diminished 
with  decreased  evidence  of  ablation.  Evidence 
tends  to  indicate  that  electrode  erosion  is 
unavoidable  for  a  starved  anode;  however,  through 
spatial  variation  of  the  mtereiectrode  spacing 
local  ohmic  heating  can  be  controlled  to  minimize 
ablation  of  the  other  components,  such  as  the 
cathode  and  backplate  insulator. 

Nomenclature 

a  magnetoacoustic  speed 
8  magnetic  induction  field 
S  electric  field 
F  thrust 
J  total  current 
J  current  density 
m  mass  flow  rate 
R  gas  constant 
rj  anode  radius 
r-  cathode  radius 
T  plasma  temperature 
u  mean  mass  velocity 
3  Hall  parameter 
v  specific  heat  ratio 
u,  permeability  of  free  space 
j  plasma  scalar  conductivity 

I nt  roduc  t ion 

The  goal  of  the  designer  of  an  MPD  arcjet  is 
to  obtain  desirable  overall  performance.  To 
achieve  this  goal,  it  is  'useful  to  manipulate 
internal  distributions  of  current  density  and  mass 
flow.  In  a  previous  paper,  Mart  Inez- Sanchez  and 
Heimerdinger  derived  an  approximate 

two-dimensional  hydrodynamic  MPD  arcjet  model 
which  indicated  that  the  performance  and  the 
internal  distributions  of  the  previously  described 

This  work  was  supported  by  the  Air  Force  Office  of 
Scientific  Research  under  grant  number 
AFOSR-86-01 19D . 


quantities  can  be  manipulated  by  an  axial 
variation  of  the  mtereiectrode  separation.'  This 
paper  presents  preliminary  lata  from  an  experiment 
based  on  the  analysis  from  reference  1. 

MPD  arcjets  can  be  distinguished  from  other 
electric  thrusters  by  their  use  of  the  jxB  Lorent,’ 
force  as  their  primary  source  of  thruat.  An 
approximate  overall  thrust  equation  shows  this 
feature  by  the  quadratic  behavior  of  the  thrust 
with  the  applied  current  and  a  weak  dependence  on 
geometry  '• 
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For  a  given  thruster  geometry  and  mass  flow  rate, 
the  efficiency  also  increases  with  current  because 
of  the  transition  from  electrothermal  acceleration 
due  to  the  ohmic  heating  of  the  plasma  to 
primarily  an  electromagnetic  acceleration. 
Unfortunately,  an  operational  limit  called  ’’onset" 
is  reached  where  both  erosion  of  the  thruster  and 
overall  voltage  have  been  noted  to  increase 
dramatically.  The  term  onset  originated  from  the 
appearance  of  megahertz  terminal  voltage 
oscillations  which  quickly  grow  in  amplitude. 
Rudolph  defined  onset  as  that  level  where  voltage 
oscillation  reach  10X  of  the  average  terminal 
voltage.1  In  a  given  geometry,  the  onset  limit 
scales  approximately  with  the  parameter  J1  m  and 
thus  correlates  with  the  exhaust  velocity.4 

In  reference  1,  the  appearance  of  onset  is 
attributed  to  an  anode  starvation  phenomenon 
caused  by  axial  deflection  of  the  local  current 
vector  by  the  Hall  effect.  In  figure  1,  the 
Lorentz  force,  which  is  orthogonal  both  to  the 
magnetic  induction  field  and  the  current  density 
vector,  has  a  component  directing  the  plasma 
towards  the  cathode.  As  the  plasma  adjacent  to 
the  anode  becomes  rarefied,  deflection  of  the 
current  vector  is  exacerbated  as  the  Hall 
parameter  grows  due  to  the  local  decrease  in 
col  1  is lonal l ty .  The  region  adjacent  to  the  anode, 
increases  in  voltage  to  attract  electrons  and 
ensure  current  flow.  Eventually  the  magnitude  of 
this  potent lsl  grows  large  enough  for  electrical 
breakdown  to  occur  which  allow*  concentrated  arcs 
to  form  on  the  anode.  This  effect  is  intrinsic  to 
MPD  arcjet  operation  unless  an  external  electric 
field  is  applied  to  halt  the  growth  of  axial 
current . 

In  addition  to  plasma  depletion  at  the  anode 
at  high  currents  and  low  mass  flow  rates,  the  MPD 
arcjet  can  operate  at  a  significant  magnetic 
Reynolds  number  which  implies  a  tendency  towards 
convection  of  the  magnetic  field  with  the  plasma.1 
Such  convection  tends  to  produce  strong 
concentrations  of  current  at  the  entrance  and 
exit,  due  to  the  presence  of  a  low  back 


electromotive  force  (emf).  Fro*  Ohs'*  law, 


]  *  JkS  -  oil  ♦  uxfl) 

the  please  suet  take  Che  entire  electric  field  in 
Che  absence  of  bade  emf,  thus  forcing  strong 
current  concent  ret  ions  to  occur  Because  of  this, 
local  heating  of  the  insulator  and  electrodes 
occurs,  accounting  for  their  associated  ablation. 

The  aia  of  this  paper  is  to  present  data  froa 
an  experiment  based  on  a  theory  that  claims  the 
ability  to  control  the  distribution  of  the 
magnetic  field,  hence  .-urrent  lenaity,  m  an  MPD 
arejet  flannel  through  areful  axial  variation  of 
the  interelectrode  separation 

The  S xper laen tal  Appar at  ua 

A  joint  experiment  between  MIT  and  RDA 
easnmgton  Research  Laboratory  m  Alexandria, 
Virginia  was  conducted  at  •he  SDA  facilities.  A 
400  xj  apacitor  bank  was  modified  into  an  eight 
stage  LC  iadder  network  to  provide  1  flat 
luasi  steady  '-urrent  pulse  to  tne  MPD  arejet.  The 
banx  has  a  total  apacitance  of  2-0  mF  where  each 
section  is  separated  by  a  10  uH  solenoid  inductor. 
The  pulse  forming  network  PFN  is  lesigned  with  a 
high  voltage  ignitron  switch  at  each  stage  so  that 
the  output  waveform  can  be  modified  depending  on 
the  sequencing  of  ignitrons  refer  to  figure  2:. 
If  a  single  ignitron  ignitron  #8  is  fired,  the 
resulting  waveform  has  a  duration  of  about  0.8  ms . 
This  system  requires  a  matched  load  of  0.2  Q 
providing  a  maximum  safe  current  pulse  of  32  kA  to 
the  MPD  arejet.  If  ignitrons  #1  and  *8  are  fired 
simultaneously,  the  impedance  of  the  bank  drops  in 
half  and  a  pulse  of  about  0.45  ms  is  achieved. 
With  this  configuration,  a  maximum  safe  current 
level  of  65  kA  can  be  supplied  to  the  arejet. 

The  PFN  is  connected  to  the  arejet.  which  is 
situated  in  a  cylindrical  stainless  steel  vacuum 
tank  6  a  long  and  0.6  m  in  diameter.  The  vacuum 
is  maintained  at  a  static  level  of  '  10' '  torr  by 
two  six  inch  diffusion  pumps.  The  tank  and 
associated  pumps  are  electrically  floating  tr 
prevent  spurious  alternate  current  paths  for  the 
PFN-arcjet  circuit.  The  MPD  arejet  assembly 
drawing  in  figure  3  shows  the  relation  of  its 
various  components.  The  exterior  of  the  arejet  is 
constructed  as  an  anode  sleeve  which  is  isolated 
from  the  vacuum  tank  by  a  Plexiglas  flange.  The 
sleeve  and  the  outer  portion  of  the  anode  are 
hardcoat  anodized  with  a  2  mil  layer  to  inhibit 
current  attachment  beyond  the  anode  lip.  The 
cathode  and  mass  injection  assemblies  both  slide 
into  the  anode  sleeve.  A  boron  nitride  insulator, 
seated  in  a  Plexiglas  assembly,  isolates  the  anode 
from  the  cathode  while  a  Mylar  film  and  tape 
jacket  surround  the  cathode  sleeve  completing  the 
anode/cathode  insulation.  By  constructing  the 
feeds  to  the  arejet  m  this  coaxial  manner,  the 
parasitic  inductance  of  the  device  is  kept  to  a 
minimum. 

Mass  is  provided  to  the  channel  from  a  large 
reservoir  placed  close  to  the  MPD  arejet  channel. 
Six  valves  feed  a  preexpansion  chamber  directly 
upstream  of  48  choked  orifices.  The  gas  flows 
through  the  choked  orifices  and  expands  through  a 
boron  nitride  insulator,  which  also  has  48  larger 
diameter  holes,  into  the  channel. 

The  anode  is  made  of  aluminum,  primarily  for 
its  machinabi 1 l ty ,  low  weight,  and  low  cost.  The 
cathodes  are  made  from  copper  primarily  for  its 
electrical  properties,  and  it  is  both  less 


expensive  and  easier  to  machine  into  the  required 
geometries  than  other  more  desirable  materials. 

A  Tektronix  data  acquisition  systaa  is 
utilited  for  data  collection,  manipulation,  and 
analysis.  It  is  kept  inside  a  large  Faraday  caga 
for  shielding  against  stray  electromagnetic  noise. 
A  block  diagram  of  the  MPD  thruster  system  is  seen 
in  figure  4. 

Initiation  of  the  diacharge  occurs  after  the 
bank  is  charged  to  the  appropriate  operating 
value.  At  this  point,  the  valves  are  opened  with 
a  3  kV  pulse  and  are  held  open  until  the  required 
steady  state  mass  flow  rate  is  reached.  After 
this  delay  approximately  10  as )  ,  the  required 
ignitrons  are  triggered,  several  kilovolts  appear 
across  the  electrode  gap,  breakdown  ensues,  and, 
after  a  short  transient  period  typically  150  us,, 
quasi-steady  operation  is  reached. 

Thruster  voltage  is  measured  by  a  1000  1 
Tektronix  voltage  probe  optically  isolated  froa 
the  data  acquisition  system.  Floating  potential 
is  measured  relative  to  the  anode  from  a  simple 
Langmuir  probe  through  a  202:1  voltage  divider. 
The  floating  potential  probe  is  constructed 
primarily  froa  a  short  piece  of  thin  tungsten 
wire.  The  total  current  is  measured  by  a  large 
Rogowski  ioop  surrounding  the  cathode  sleeve  at 
the  base  of  the  power  connection.  This  probe 
measures  the  time  rate  of  change  of  the  enclosed 
magnetic  flux,  which  is  integrated  to  give  the 
magnetic  field,  and,  froa  Ampere’s  law,  the  total 
enclosed  current.  An  active  integrator  with  a 
377  us  integration  time  constant  and  a  70  ms 
associated  droop  time  is  used  with  this  Rogowski 
loop.  Enclosed  current,  at  points  inside  the 
plasma,  is  measured  froa  a  small  10  turn  Rogowski 
loop  enclosed  in  a  long  F*yrex  tube.  Again,  an 
active  integrator  with  a  3.49  ms  integration  time 
constant  and  a  70  as  associated  droop  time  is  used 
to  integrate  the  time  derivative  of  the  magnetic 
induction  field. 

For  this  experiment,  the  bank  was  fired  with 
#1  and  *8  ignitrons  so  that  an  operating  point  of 
60  kA  at  4.0  g/s  of  Argon  was  obtained. 

Two  MPD  arejet  channels  were  constructed  and 
diagnosed.  A  constant  area  channel  consists  of  a 
cylindrical  anode  and  concentric  cathode  with  a 
constant  interelectrode  separation  of  1.8  cm.  A 
second  channel,  the  modified  flared  channel,  was 
based  on  a  calculation  froa  the  theory  in 
reference  1  for  uniform  current  density 
distributions  along  the  electrodes.  To  minimize 
cylindrical  effects,  the  annular  channel  geometry 
was  based  on  a  large  cathode  diameter  compared  to 
the  mterelectrode  gap.  For  a  mass  flow  rate  of 
4.0  g/s  of  Argon  at  an  operating  point  of  about 
60  kA,  a  channel  with  the  area  variation  shown  in 
figure  5  is  predicted  by  the  theory.  In  order  to 
examine  the  behavior  of  the  two  channels,  we  chose 
to  operate  in  the  vicinity  of  onset,  where  onset 
is  predicted  by  the  theory  as  the  solution  where 
zero  density  is  encountered  somewhere  along  the 
anode.  The  design  theory,  for  which  this 
calculation  is  based,  neglects  the  initial 
ionization  processes  that  occur  in  a  region  of  the 
discharge  dominated  by  ohmic  heating  where  large 
changes  in  plasma  velocity  and  density  occur.' 
For  this  reason,  a  low  speed  near  the  entrance  for 
a  given  mass  flow  rate  leads  to  large  variations 
in  the  mterelectrode  separation.  This  is 
iMediately  seen  in  figure  5  at  the  channel 
entrance  upstream  from  the  point  of  minimum  area. 
However,  for  this  experiment,  in  order  to  ensure 
reliable  ignition  at  the  inlet,  only  the  portion 


downstream  froa  th#  point  of  uniaia  area  wm a 
flared.  The  upatreaa  portion  remains  fixed  at 
1.8  cm.  Hence  thia  device  la  called  the  modified 
flared  charnel.  In  both  geometries,  the  anode 
remain#  the  same  and  ail  channel  variation  1a 
taken  by  the  cathode.  This  la  done  for  simplicity 
in  the  machining  of  the  electrodes. 

Results  and  Discussion 

Both  channels  were  probed  to  determine 
current  and  floating  potential  distributions. 
Typical  probe  responses  are  seen  in  figures  6  and 

The  enclosed  current  and  floating  potential 
contours  for  the  two  cnannels  are  found  in  figures 
8  and  9  respectively.  Determination  of  the  plasma 
potential  froa  floating  potential  measurements 
necessitates  knowledge  of  the  electron  temperature 
and  and  ionic  species."  Other  work  has  shown  that 
this  temperature  variation  is  small  so  that  the 
J 1  iterances  in  floating  potential  is  sufficient 
tor  the  calculation  of  electric  fields. 

In  both  channels,  large  axial  currents  are 
r  'und  in  tne  anode  region.  This  implies  large 
.oral  va.ues  of  the  Hall  parasieter.  From  the 
t.oating  potential  contours,  there  aiso  appears  to 
Se  a  large  positive  potential  drop  at  both  anodes. 

A  large  Hail  parameter  implies  a  a  deficiency  of 
■large  arriers  which  orreiates  with  the  large 
positive  anode  fall  required  to  pass  all  of  tne 
current.  Based  on  reference  1,  it  is  apparent 
'hat  Poth  cnannels  are  operating  Beyond  the  onset 
level . 

Typical  terminal  voltage  traces  do  not 
tispiay  the  high  frequency  oscillations  that 
typically  characterize  onset.  Megahertz  voltage 
oscillations  are  seen,  however,  in  the  floating 
potential  charact er 1st lcs .  Terminal  voltage 

scaling  shown  in  figure  10  does  not  display  the 
cubic  dependence  of  voltage  that  one  would  expect 
to  precede  onset.1  * 

Turchi  suggests  that  an  increase  in  ablation 
as  the  exhaust  speed  exceeds  the  Alfven  critical 
speed,  based  on  ablation  energy,  may  preclude  a 
unic  dependence  of  voltage  on  current  in  MPD 
arcjet  operation,  even  if  the  dominant 
acceleration  mechanism  is  electromagnetic.’ 

Examination  of  both  thrusters  shows  strong 
evidence  of  ablation.  Both  anodes  show  erosion 
primarily  at  their  exit  lip.  Both  cathodes  also 
show  erosion,  especially  in  the  entrance  region, 
in  fact,  after  less  than  100  discharges,  copper 
from  the  cathode  is  plated  on  the  boron  nitride 
insulating  injector  plate.  The  anode  wear  is 

similar  m  both  geometries,  however,  the  cathode 
wear  is  markedly  different.  For  the  constant  area 
channel,  bands  about  a  half  a  centimeter  wide  are 
seen  at  both  the  entrance  and  the  exit.  The  band 
at  the  entrance  looks  like  freshly  machined 
upper,  indicative  of  strong  and  even  erosion. 
The  band  at  the  exit  is  well  defined  but  is  not 
juite  as  clean  as  the  one  at  the  entrance.  The 
band  at  the  entrance  are  largely  composed  of 
semicircular  clean  regions  adjacent  to  the  mass 
.ejection  sites.  More  clean  spots  are  noted  about 
a  half  to  one  centimeter  further  downstream 

iirectly  in  line  with  the  mass  injection  holes  in 
■he  boron  nitride  insulator.  The  remaining 
athode  shows  strong  evidence  of  fine  line  arc 
structures  extending  from  the  cathode  downstream 
“nd  almost  to  the  cathode  root.  These  fine 

structures  are  branch-1  ike  and  tend  to  branch  in 
an  upstream  direction. 

In  the  modified  flared  channel,  the  clean 


band  at  the  exit  la  absent ,  but  tbsrm  is  a  widm 
clean  band  at  the  entrance  between  0.75  cm  to 
1.3  cm  in  extent.  Unlike  the  other  band*,  this 
one  has  a  very  poorly  defined  downstream  edge. 
The  fine  line  structure#  are  completely  absent  on 
this  cathode.  Drawings  of  both  cathodes  are  found 
in  figures  11  and  12. 

Figure  13  shows  a  graph  of  the  current 
density  along  both  cathode  surfaces.  In  both 
cases,  strong  concentrations  are  found  in  the 
initial  portion  of  the  channel.  In  the  constant 
area  channel,  the  current  density  drops  to  a  low 
value  at  aid-channel  and  then  increaaes 
dramatically  near  the  exit.  The  current  density 
in  the  modified  flared  channel  remains  high  until 
just  before  the  beginning  of  the  flared  region 
where  it  subsequently  drops  significantly  below 
the  level  found  in  the  constant  area  channel. 

The  theory  tends  to  support  the  trends  seen 
in  figure  13.  At  the  channel  entrance,  the  plasma 
is  moving  slowly  so  the  back  emf  is  low.  In  both 
channels,  the  interelectrode  distance  at  the 
entrance  is  the  same,  and  both  plasmas  must  taxe 
the  total  electric  field  in  the  absence  of  any 
significant  back  emf.  At  the  exits  the  magnetic 
field  drops  towards  zero  and  once  again  the  back 
emf  is  small.  In  this  region,  the  plasma  once 
again  must  take  the  total  electric  field;  however, 
in  the  modified  flared  channel,  the  local  electric 
tield  decreases  due  to  the  increased 
interelectrode  separation,  so  a  strong  current 
concentration  is  mitigated.  Since  the  current 
density  is  lower  in  the  flared  region,  and  the 
same  current  is  passed  as  in  the  constant  area 
channel,  the  current  density  in  the  upstream 
portion  of  the  modified  flared  channel  must 
increase.  This  may  account  for  higher  ohmic 
heating  which  may  lead  to  an  elevated  plasma 
temperature,  and  ther"fore  higher  ionization 
fraction  and  electrics.  onductivity,  in  the 
modified  flared  channel. 

The  theory  also  indicates  that  high  magnetic 
Reynolds  number  flows  in  constant  area  channels 
tend  to  have  strong  current  concentrations  at  the 
exit  as  the  magnetic  field  goes  to  zero.  Also  it 
shows  that  the  magnetic  field  in  a  flared  channel 
is  more  evenly  distributed  due  to  the  controlled 
plasma  expansion  and  the  convection  of  the 
magnetic  field  with  the  fluid  at  higher  magnetic 
Reynolds  numbers.  In  this  operating  regime,  the 
magnetized  plasma  behaves  analogously  to  an 
ordinary  compressible  gas  except  that  the  speed  of 
sound  is  replaced  by  the  magnetoacoustic  speed, 
defined  as 

B1  1  ■  1 

a  =  t tRT  ♦  - — ) 

which  is  a  combination  of  the  speed  of  sound  and 
the  Alfven  speed. 

Data  for  the  constant  area  channel  does  not 
seem  to  exhibit  this  behavior.  The  theory 
neglects  the  effect  of  viscosity  which  may  account 
for  some  of  the  differences.  In  this  channel,  the 
presence  of  a  large  axial  component  of  the  current 
acts  to  deplete  the  anode  and  pressurizes  the 
cathode.  The  moving  plasma,  blocked  by  the 
cathode  must  move  along  it,  which  may  be 
sufficient  for  viscous  dissipation  to  slow  the 
plasma.  If  the  plasma  speed  is  sufficiently 
lowered,  the  corresponding  magnetic  Reynolds 
number  for  the  arcjet  may  be  low  enough  for  a  more 
uniform  axial  decrease  of  magnetic  field.  If  the 
magnetic  field  is  sufficiently  low  at  the  exit, 
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ordinary  I u  dynaaic  choking  can  occur  os  in  tha 
caaa  of  ordinary  gaa  dynaaica  with  Tanno  flow  for 
flow*  with  fnctioo  and  Rayleigh  flow  for  flow* 
with  heating),  thus  lieiting  the  effective  exhaust 
velocity  to  that  corresponding  to  Mech  1. 

However,  the  aodified  flared  channel  does 
exhibit  ease  of  the  trends  anticipated  by  the 
theory.  Although  preeaunzat  ion  of  the  cathode  is 
evidenced  by  the  current  contours,  the  flaring  of 
the  cathode  tends  to  reduce  the  obstruction  of  the 
plaaaa  action.  Also,  an  increased  plaaaa 
teaperature  can  lead  to  a  strong  drop  in  plaaaa 
viacoeity  as  shown  in  figure  14. • 5  This  could  act 
to  sufficiently  lower  the  viscous  losses  expected 
in  this  region.  Therefore,  the  sagnetic  Reynolds 
nuaber  say  not  be  reduced  significantly,  in  fact, 
the  thermodynaaic  expansion  encountered  in  the 
nozzle  region  would  allow  for  an  increase  in  the 
plaaaa  velocity  beyond  the  appropriate  sonic  or 
magnetoacoust ic  speed,  which  would  translate  to  an 
increase  in  the  aagnetic  Reynolds  nuaber. 

This  is  further  evidenced  by  the  difference 
m  the  terainal  voltages  measured  in  the  arcjets. 
Mart inez-Sanchez  shows  that  the  terainal  voltage 
decreases  with  increasing  aagnetic  Reynolds 
nuaber.*  In  this  experiaent,  the  constant  area 
cnannel  is  seen  to  operate  at  a  Bean  voltage  of 
'8.1  t  1.6  volts,  while  the  modified  flared 
:hannel  operates  at  70.0  -  1.8  volts. 

The  plaaaa  discharge  near  the  anode  appears 
quite  different  froa  the  cathode.  The  anode 
current  densities  and  voltage  drops  are  plotted  in 
figure  15.  Both  geoaetries  exhibit  large 
potential  drops  indicative  of  starved  anodes.  for 
the  aost  part,  the  magnitude  of  the  voltage  drop 
follows  the  aagnitude  of  the  current  density.  The 
only  place  where  there  is  a  deviation  froa  this 
behavior  is  at  the  entrance  to  the  aodified  flared 
channel.  Here,  high  current  density  la  associated 
with  a  low  voltage  drop.  This  tends  to  suggest 
the  development  of  concentrated  arcs, 
substantiated  through  examination  of  the  anode. 

In  this  initial  region,  there  is  a  high 
Jeaand  to  pass  current  due  to  the  low  back  eaf, 
but  the  gas  has  a  low  conductivity  since  it  has 
not  yet  become  fully  ionized.  The  local  electric 
fields  may  not  yet  be  sufficient  to  support  a 
diffuse  discharge  so  the  plasma  may  have  no  other 
choice  but  to  arc  strongly  and  locally  heat  the 
plasma.  As  the  plasma  is  convected  downstreaa, 
its  conductivity  increases,  and  the  arc  may  become 
more  diffuse.  Eventually,  no  severe  arcing  is 
necessary,  as  the  conductivity  and  the  back  eaf 
nave  risen  to  a  level  so  that  a  more  diffuse 
lischarge  may  be  supported.  This  is  evidenced  by 
the  wear  pattern  on  the  inside  of  the  anode. 
Arcing,  adjacent  to  the  injectors  at  the  anode,  is 
seen  with  tracks  that  appear  to  continue 
downstreaa.  As  the  arcs  continue  downstreaa,  they 
appear  to  broaden  and  weaken  until  the  wear  on  the 
anode  surface  becomes  more  azimuthal ly  uniform; 
►his  occurs  about  a  third  of  the  way  down  the 
anode  froa  the  channel  entrance. 

Significant  erosion  can  occur  in  other  areas 
of  the  AW’D  arcjet  for  reasons  that  are  not 
necessarily  associated  with  anode  current 
conduction  onset.  In  this  experiaent,  erosion  of 
the  cathode  has  been  locally  decreased  by  flaring 
'he  cathode.  Theory  suggests  that  this  is  also 
valid  for  the  entrance  region.1  By  decreasing  the 
local  electric  field,  the  current  is  spread  out, 
thereby  decreasing  its  density  and  the  associated 
local  ohaic  heating.  This  not  only  protects  the 
cathode,  but  the  insulator  region  as  well. 


Conclusion 

Preliminary  examination  of  data  froa  an  MPO 
arcjet  experiaent  has  shown  that  arcjet  operation 
can  be  modified  through  a  variation  of  the 
interelectrode  separation.  Further 

experimentation  should  be  done  to  study  the  effect 
of  increasing  the  interelectrode  distance  at  'he 
entrance  in  a  siailar  Banner  as  was  done  at  the 
exit  to  alleviate  the  high  concentration  of 
current  along  the  cathode.  If  the  expected 

decrease  in  current  iensity  is  encountered. 
Additional  experimentation  using  better  cathode 
materials  with  superior  ablation  properties  suih 
is  thoriated  tungsten  or  barium  oxide  impregnated 
tungsten  ■  :  should  be  lone,  as  well  is  ont  inued 
study  of  the  MPD  arc  ;et  plasma  properties  so  is  'o 
understand  loss  mechanisms  such  as  friction  ,n  'he 
MPD  arcjet  And  finally,  idditional 

experimentation  should  be  conducted  to  study  'be 
'nruster  operation  oeiow  onset. 
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To  runiK  fxmibk  nkaot  velocity  li™  Hallow  ia  high  power  mgartotiasaadyaaink  (MPD)  arejets.  ■ 
power  balioce  aiodei  is  coutrwcled.  A  rpecnU  fee  Is  re  of  this  model  is  tkai  Ike  dnckirye  Ikkkneu  is  Iwveraeiy 
pro  port  toul  to  (ke  velocity  of  Ike  kigk-speed  elect  ncsBy  CO  ad  acting  flow  crested  ky  Ike  dtsefcarge  sod  its  self 
magnetic  field.  Resistive  dissapstios  Ikes  sales  wttk  fleet  ro  mx I  set  k  tknasa  swd  flow  speed.  AMaiioa  of  arc  yet 
materials  adds  to  Ike  total  mass  flow  wkes  Ike  iapal  propeilaat  flow  is  iaxafficieat  to  carry  away  Ike 
dissipated  power.  At  ktfk  power  levels.  Ike  exkaast  velocity  may  pialeaa  keaase  of  sack  mass  addltioa  at 
values  comparable  to  Alfvea  critical  speed.  At  lower  powers,  iacompiete  ioaizatioa  aad  coupllag  of  Ike  ia- 
jected  aeairal  now  in  Ike  discharge  as  also  result  ia  sack  a  velocity  piatraa  fas  taterpreted  by  votta(e  aad 
pluaie  measurements!  The  vottage-carmi  ckaract ensues  of  MPD  arejets  are  discussed  aad  rears rrk  direc¬ 
tions  are  tad  tea  ted 


Nomenclature 

=  magnetic  (induction)  field 
=  first  ionization  level  of  propellant  atom 
=  geometric  factor  relating  discharge  volume  to 
thruster  area 

-  current  density 
=  total  current 

=  total  current  for  which  3  =  1 

-  mass  of  propellant  atom 
=  total  mass  flow  rate 

-  mass  flow  rate  of  ablated  material  and  propel¬ 
lant  material,  respectively 

=  heat  absorbed  per  unit  mass  of  ablated 
material 

-  heat  absorbed  per  unit  mass  of  propellant 
material 

=  characteristic  radius  of  thrust  chamber 
=  ratio  of  outer-to-mner  current  attachment 
radii 

=  magneto  Revnolds  number 

-  exhaus'  speed 

-  Alfven  critical  speed 
characteristic  speed  ( for  d  -  I) 

-exhaust  speed  of  ablation  thruster  im,m0) 
speed  of  electrically  conducting  ponton  of  flow 
-electron  flow  velocity 

ideal  exhaust  speed  (full  coupling,  no  ablation) 

-  voltage 

-  mass  density 

-  electrical  conductivity 
=  magnetic  permeability 

-  characteristic  discharge  thickness 

=  characteristic  discharge  thickness  for  ablatton  arc 

-  characteristic  volume  of  discharge 

-  ratio  of  ablated  mass  flow  rate  to  input  pro¬ 
pellant  mass  flow  rate 

-  ablation  parameter  indicating  relative  ease  of 
ablating  material.  =  Q ,  >Q,  (high  o  means  abla 
tion  is  difficult) 
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3  =  operating  parameter  indicating  ideal  thrust  power 

relative  to  convective  capability  of  input  mass 
flow,  =  (ur'ua):  (high  8  means  significant  excess 
thermal  power) 

r  =  fraction  of  excess  thermal  power  (dissipated 

power  minus  power  convened  by  propellant 
flow)  causing  ablation 

Introduction 

FOR  several  years,  it  has  been  suggested  that  the  exhaust 
velocity  of  a  magnetoplasmadynamic  (MPD)  arejet  is 
limited  by  the  so-called  Alfven  critical  speed1: 

ken,  =  (2£/m)  ’’  (1) 

where  E  is  the  first  ionization  level  of  a  fuel  atom  of  mass 
m  A  few  theories,  based  on  quite  different  physical 
models,'4  have  led  to  this  speed  as  a  characteristic  value  of 
the  exhaust  velocity  at  which  electrode  erosion,  insulator 
ablation,  and'or  "plasma  instabilities"  develop  that  serve  lo 
limn  arejet  performance  In  particular,  for  arejet  currenl  J 
and  mass  flow  m,  the  "onset"  of  difficulties  appears  to  oc¬ 
cur  ai  values  of  P/m  that  result  in  exhaust  speeds  near 
Alfven  critical  speed.  The  present  discussion  invokes  the 
electromagnetic  structure  of  the  discharge  flow  and  a 
requirement  for  power  balance  in  steady  state  in  order  to 
derive  self-consistent  mass  flow  rates  and  exhaust  speeds. 
Implications  for  voltage-current  characteristics  are  described 
and  further  research  directions  are  suggested. 

Model 

The  basic  equation  for  the  current  distribution  in  a 
longitudinal  plasma  flow  is 


D,  /  B  \ 

Dr  -o'  nan 


where  p  is  mass  density,  o  electrical  conductivity  (assumed 
constant),  and  convective  derivative  is  based  on  the  electron 
fluid  velocity  u,  (For  high-density  flows,  the  plasma  flow 
velocity  u  may  be  substituted  for  u,  )  By  dimensional 
analysis,  the  characteristic  scale  over  which  the  magnetic 
field  varies  ihrough  the  plasma  flow  in  steady  state  is 

6  =  I  /<t|iu  (3) 
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Without  Ho*,  the  magnetic  field  would  decrease  linearly  be¬ 
tween  the  upstream  and  downstream  boundaries  of  the 
thruster  As  the  conducting  flow  accelerates,  it  convects 
magnetic  flux  away  from  the  upstream  region  and  concen¬ 
trates  llus  near  the  downstream  (fl  =  0)  boundary  A  bifur¬ 
cated  current  density  distribution  is  thereby  obtained  Higher 
flow  speeds  result  in  greater  current  concentration,  and  more 
intense  dissipation  The  consequence  of  such  dissipation  is 
increased  ablation  of  electrodes  and  insulators  Ablation, 
however,  contributes  to  the  total  mass  flow  rate  and  thereby 
lowers  the  flow  speed,  reducing  the  current  concentration.  A 
self-consistent  ,olution  should  therefore  be  expected. 

As  a  simple  model,  let  convection  provide  the  primary 
mechanism  balancing  resistive  dissipation,  with  ablation  ab¬ 
sorbing  the  heat  not  carried  away  by  the  fuel  mass  flow 

'■”<(?«  =>l  ~mFQf]  (4) 

where  „■  is  the  current  density,  v  the  volume  in  which  dissi¬ 
pation  occurs,  m  the  mass  flow  rale,  and  Q  the  heat  ab¬ 
sorbed  per  unit  mass  Subscripts  ,4  and  F  refer  to  ablated 
and  tuel  mass,  respectively  The  factor  t  allows  energy  to  be 
lost  without  causing  ablation  leg.,  radiation  to  free  space, 
thermal  conduction). 

The  current  density  may  be  written  in  terms  of  the  total 
current  J  and  the  characteristic  scale  a,  which  also  provides 
the  volume,  so 

L/‘*/t»)  =  (l  aHJ  2rrd):Tr~6g  =  (tl/4x)J2ug  (5) 

where  ?  is  a  geometric  factor  that  depends  on  the  particular 
thruster,  and  u  the  exhaust  speed.  The  electromagnetic  thrust 
provides  this  speed 

u-[n;4iuf-  m  )|  £*(/-,  r;  )-*>'.  |  (6) 

where  m  is  the  total  mass  flow  rate  mf  ♦  m4.  Substitution  in 
the  equation  for  ablation  results  in  a  quadratic  equation: 

im,  -  m,  i  \m,Q,  -  mt  (£>,  « )  ] 

-  1  w  4r  i-  J*g  |  ( fa  ( r.  r, )  --  )  j  (7) 

The  ablation  rate  relative  to  the  fuel  flow  rate  is  then 
£T=  1  m,  )  -  !  !  to  *  l);  -  4o(d  -  D]  -  (a*  1)  I  la  t8) 

where  a  Q ,  >Q,  and 

d  =  '  I  a  4r  )  t  p  mf  )  I ;  [g (  6cT.  r  -  *4  )]  Qf  (9) 

For  del.  -here  is  no  physical  solution  for  the  model  as 
constructed  There  is  insufficient  dissipation  to  account  for 
the  heat  that  could  be  carried  away  by  the  injected  mass  flow 
rate  absorbing  an  energy  per  unit  mass  QF  This  insuffi¬ 
ciency  can  be  resolved  simply  by  reducing  Q, ,  i.e .,  allowing 
incomplete  ionization  of  the  injected  flow.  For  d>  I,  the  fuel 
mass  flow  cannot  cool  the  thruster  sufficiently  and  ablation 
occurs  The  degree  of  ablation  and  thus  the  exhaust  velocity 
attained  for  a  given  current  flow  depends  on  the  heat  that 
.an  be  absorbed  by  the  ablated  material,  expressed  relative 
to  the  fuel  by  the  parameter  a. 

The  condition  J  -  I  provides  a  value  for  p  m,  for  which 
m  a  -  0  and  the  dissipation  is  exactly  balanced  by  convection 
in  the  fuel  mass  flow  This  value  of  P  m,  substituted  in  the 
equation  for  exhaust  velocity  gives 

u  -  u,,  =  IQ,g<l*s:  r,  -  ' . )  |  (10) 


for  Qf-E  m  and  g<6vr;,rl  ♦  *4 )«2,  Thus,  cor¬ 

relation  of  self-field  plasma  thruster  performance  with 
Alfven  critical  speed  might  be  expected  simply  on  the  basis 
of  ablation  when  the  fuel  flow  rate  is  insufficient  for  thruster 
cooling 

In  Fig.  I.  the  actual  speed  (i.e..  with  ablation)  is  compared 
with  u,.  the  ideal  exhaust  speed  tno  ablation),  as  the  ideal 
speed  is  increased  relative  to  u0  If  there  is  no  possibility  of 
ablation  (<j=od),  then  the  exhaust  speed  can  increase  in¬ 
definitely.  In  the  presence  of  plasma  at  1-3  eV,  however, 
solid  surfaces  can  be  expected  to  ablate  and  ionize  about  as 
readily  as  neutral  fuel  gas  (a»l).  Thus,  a  little  excess 
heating  is  rapidly  compensated  by  additional  mass  flow, 
resulting  in  an  exhaust  speed  “plateau"  at  about  u0 »  u^, 

For  a  =  0.  which  would  correspond  to  an  ablation-fed 
thruster.  Fig  I  indicates  that  the  exhaust  sjyeed  decreases 
relative  to  u0  based  on  Qr.  Proper  substitution  of  this 
limiting  case  (mr  =  0)  in  the  ablation  cooling  model,  how¬ 
ever,  provides  a  constant  operating  speed 

“4  =  f <2w« ( £»< ^2/ M  |  1  HD 

that  is  basically  equivalent  to  operation  at  d  =  1  The 
characteristic  discharge  thickness  for  an  ablation  arc  is  then 

6,  =  ( 1/oiiu.,  )  (12) 

Incomplete  Ionization  and  Coupling 

As  noted  earlier,  the  situation  of  I  can  be  associated 
with  incomplete  ionization  of  the  input  mass  flow  If  the 
mean  free  path  for  momentum  exchange  between  ions  and 
neutrals  is  small  compared  to  the  discharge  and  thruster 
dimensions  and  the  input  flow  is  well  mixed  into  the  dis¬ 
charge,  then  it  would  be  appropriate  to  divide  the  electro¬ 
magnetic  thrust  by  the  full  mass  input  rate.  The  actual  ex¬ 
haust  speed  should  then  follow  the  ideal  speed  as  indicated 
in  Fig.  I  for  0<l.  With  poor  mixing,  however,  or  at  low 
densities,  the  10ns  can  slip  relative  to  the  neutrals,  and  the 
velocity  of  the  electrically  conducting  fluid  can  be  higher 
than  the  ideal  speed.  In  the  limit  of  complete  slip,  the  mass 
flow  subjected  to  the  electromagnetic  thrust  scales  with 
dissipation.  The  situation  is  completely  analogous  to  the 
pure  ablation-fed  thruster,  with  electrical  dissipation  effec¬ 
tively  'ablating''  the  input  neutral  mass  flow  The  conductor 
speed  will  be  given  by: 

uc  =  I  Q,g[t.(r:>r,  )  4-  >4  ]  |  '  (13) 


Fig.  I  Viriattoa  of  actual  speed  with  expected  speed,  both  relative 
to  critical  speed  u,.  as  a  faaciloa  of  energy  pee  aait  at  ass  foe  able 
lloa  relative  to  fact,  a  ■  Q4  tQt  (labels  refer  to  reflates  of  opera- 
fioa  discussed  la  text.) 
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In  (he  ion  slip  limit,  the  ionized  fraction  of  the  flow  would 
achieve  a  high  speed  at  modest  current  levels  and  would  ex¬ 
hibit  constant  speed  until  d  =  1  [which  corresponds  to  com¬ 
plete  tsinglet  ionization  of  the  input  mass  flow] 

Volmge-CurTenl  Characteristics 

To  (he  extent  that  the  voltage  across  an  MPD  arejet 
largels  reflects  the  electromotive  lux  j]  effects  of  plasma 
'Tow  at  significant  magnetic  Reynolds  number,  the  arejet 
voltage-current  characteristic  should  display  regimes  of 
ipe-ation  associated  with  the  veloctt>  of  the  conducting 
(Tow  In  Fig  2,  the  components  of  the  arejet  voltage  are 
sketched  Anode  and  cathode  falls  provide  a  significant  bias 
tha:  mas  be  nearly  constant  with  respect  to  total  current 
(although  interesting  variations,  particularly  of  the  anode 
tali  might  be  expected)  In  idealized  MPD  operation,  for 
which  ihe  mass  flow  is  only  the  input  mass  flow  (no  abla¬ 
tion:  and  a  single  heavy-particle  velocity  exists  (full  toniza- 
■ion  and  or  no  ion  slip),  the  conducting  flow  speed  scales  as 
J'  so  the  voltage  l  -  uB  ~  J'  At  high  current  levels  (ff>l) 
and  limte  a,  however,  the  mass  flow  includes  ablated 
materia!  and  the  speed  increases  less  rapidly  with  current, 
pe-haps  even  becoming  limited  to  a  constant  value  The  back 
t MF  component  of  voltage  will  then  scale  linearly  with 
current 

If  tne  input  mass  flow  is  not  fulls  ionized  and  colhsion- 
alirv  between  ions  and  neutrals  is  not  sufficient  to  prevent 
significant  ion  slip,  then  the  flow  speed  that  determines  the 
ha-k  electromotive  force  (EMF)  will  be  higher  than  that 
computed  from  the  electromagnetic  thrust  divided  by  the  full 
input  mass  flow  rate  Since  this  flow  speed  will  also  be  con¬ 
stant,  scaling  with  Qf ,  the  back  EMF  voltage  will  increase 
linearis  with  J  at  a  level  above  that  for  the  idealized  f  varia¬ 
tion  Intermediate  degrees  of  collisionaltty  between  full  slip 
and  complete  ion-neutral  coupling  provide  voltage-current 
characteristics  between  the  linear  and  cubic  variations. 

If  variations  in  conductivity  and  current  path  length  with 
total  current  level  (and  mass  flow  rate)  are  ignored  compared 
to  the  basic  variation  of  discharge  thickness  with  arejet 
operation,  then  two  cases  can  be  identified  for  the  resistive 
component  of  the  voltage  For  high  magnetic  Reynolds 
number  flow,  the  discharge  thickness  scales  as  (ojzul  ,  so 
'he  resistance  will  increase  as  u  and  the  resistive  voltage  will 
.arv  av  uj  This  variation  will  therefore  be  included  in  the 
-amt  manner  as  the  back  EMF 

Ir,  some  regions  ot  the  flow,  particulars  near  the  pro- 
peltai’’  niet.  where  is  low.  the  current  distribution  mas 
cc  Je’ermined  bv  processes  such  as  ionization  and  ihermai 


*  s  9  1  i  z  1  *  1  •  1  e 


•of  2  '<kflrh  of  MPD  trcjfl  volitgr  -current  characteristics 
Voltage  *  has  arbitrary  units  and  current  J  is  relator  to  critical  cur 
'em  tor  which  u,  =  u#  ij  i*  -  u,  u9  =  £f 1  (labels  refer  to 
regimes  of  operation  discussed  in  text  i 


transport  (t.e,,  conductivity  variations  in  the  direction  of  the 
flow  are  important)  If  the  discharge  thickness  is  constant, 
for  example,  then  the  resistive  component  of  voltage  will  in 
crease  linearly  with  current  and  might  combine  with  the  back 
EMF  voltage  for  constant  flow  speed  (due  either  to  in 
complete  ionization  and  coupling,  or  significant  ablation) 

The  lota)  voltage-current  characteristic  for  an  MPD  arejet 
can  display  a  considerable  variation  depending  on  the  par¬ 
ticulars  of  geometry  (e  g.,  inlet  port  location  and  shape)  and 
material  selection  (a  value,  collision  cross-sections)  A  linear 
increase  of  voltage  with  current,  above  constant  fall 
voltages,  above  constant  fall  voltages,  may  be  expected  at 
low  magnetic  Reynolds  number  and  for  incomplete  lontza 
non  and  coupling  Transition  to  cubic  variation  of  voltage 
with  current  should  then  occur  as  full  ionization  -  1)  is  at¬ 
tained,  but  may  be  achieved  at  lower  values  of  current  if  the 
magnetic  Reynolds  number  and  degree  of  ion-neutral  colli 
sionaltty  are  both  high  This  transition,  however,  mav  be 
obscured  if  ablation  becomes  important  and  the  voltage 
returns  to  a  linear  variation  with  current  Two  quite  distinct 
aspens  of  arejet  design  (inlet  geometry  and  insulator 
material)  could  thus  combine  to  provide  a  rather  linear 
voltage-current  characteristic  for  the  range  ol  experimental 
operation 

To  the  extent  that  the  arejet  current  is  maintained  bv  the 
relatively  high  impendance  of  the  power  source,  temporal 
variations  in  the  back  EMF  component  of  the  voltage  reflect 
similar  variations  in  flow  speed  Such  flow  speed  variation 
affects  the  current  density  distribution,  i.e..  the  discharge 
thickness  (i-u  '),  At  high  power  levels,  significant  varta 
ttons  in  ablation  may  then  occur  which  in  turn  result  in  flow 
speed  variations  As  a  simple  example,  suppose  that  the 
coupling  of  dissipation  to  ablatable  surfaces  (as  expressed  by 
the  efficiency  «)  changes  abruptly  from  t  =  I  for  a  high-speed 
flow  (no  ablation)  to  <  =  0,  if  ablation  lowers  the  flow-speed 
and  thereby  broadens  the  discharge.  The  arejet  flow  field 
could  then  flicker  in  times  scaled  by  the  chamber  length 
divided  by  the  flow  speed.  Voltage  oscillations  at  a  few 
hundred  kilohertz  would  then  be  observed  (with  amplitudes 
corresponding  to  the  voltage  differences  between  different  a 
values  in  Fig.  2).  Such  oscillation  in  arejet  voltage  and 
flowfield  could  significantly  increase  the  (ranvport  of  heal 
and  mass,  further  degrading  performance 

Conclusions 

The  modeling  based  op  power  balance  in  magneto- 
plasmadynamic  (MPD)  arejet  operation  is  essennallv  a 
dimensional  argument  that  indicates  scaling  relationships  To 
predict  from  first  principles  the  detailed  behavior  of  an  arc 
jet  of  particular  geometry,  material  choice,  and  propellani 
type  would  require  a  substantially  more  complex  calcula 
tional  tool,  namely  a  two-dimensional  computer  code  The 
present  state  of  the  art  for  such  codes  is  quite  sufficient  for 
hvdrodvnamic  calculations  including  radiauon  transport, 
plasma  chemislrv.  and  interactions  with  solid  boundaries 
C  olhsionless  motion  of  ensembles  of  charged  particles  can 
also  be  computed  in  self-consistent  electric  and  magnetic 
fields  Estimates  of  collisions!  mean  tree  paths  n  arcjei 
dimensions,'  however,  indicate  that  at  least  portions  of  the 
MPD  arejet  operate  in  a  transition  flow  regime,  possible 
with  nonequilibrium  <i  e  ,  evolving)  chemislrv  Phenomena 
such  as  ion  slip  and  incompleie  ionization  and  or  mixing  ol 
ihe  injected  neuiral  flow  mav  be  inaccessible  n>  both  hvdro 
dvnamn  and  collivionlcw  codes 

In  lieu  ol  detailed  calculations,  the  present  imvdeling  ai 
tempts  to  provide  a  structure  to  support  understanding  of  ex 
perimcmal  atcjei  behavior  Experimental  measurements  of 
ion  vs  neuiral  velocities  as  a  function  ol  position  wnhm  the 
arejet  chamber  would  be  useful  to  assess  the  degree  of  colli 
sionalilv  land  incidentally  the  degree  of  lom/ationl  C  otrela 
lion  ol  collisionalits  and  ionization  with  overall  arciel 
behavior  (e  g  .  relationship  to  onsel  ol  noise  on  soilage 
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signals)  and  with  variations  in  chamber  geometry,  especially 
the  fuel  injection  arrangement,  could  then  lead  to  improved 
criteria  for  MPD  arcjet  design. 
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